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Abstract 

A relativistic hadronic model for nuclear matter and finite nuclei, which 
incorporates nonlinear chiral symmetry and broken scale invariance, is pre- 
sented and applied at the one-baryon-loop level to finite nuclei. The model 
contains an effective light scalar field that is responsible for the mid-range 
nucleon-nucleon attraction and which has anomalous scaling behavior. One- 
loop vacuum contributions in this background scalar field at finite density are 
constrained by low-energy theorems that reflect the broken scale invariance of 
quantum chromodynamics. A mean- field energy functional for nuclear mat- 
ter and nuclei is derived that contains small powers of the fields and their 
derivatives, and the validity of this truncation is discussed. Good fits to the 
bulk properties of finite nuclei and single-particle spectra are obtained. 



I. INTRODUCTION 



Descriptions of nuclear matter and finite nuclei, which are ultimately governed by the 
physics of low-energy quantum chromodynamics (QCD), are efficiently formulated using low- 
energy degrees of freedom — the hadrons. In the absence of direct derivations from QCD, 
such effective descriptions should be constrained by the underlying symmetries of QCD, 
both broken and unbroken. Nevertheless, the appropriate realization of these symmetries for 
phenomenological models is not yet established. In this paper, we explore some consequences 
of applying QCD symmetry constraints to a relativistic model of finite nuclei that features 
a light scalar meson. 

At present, the most developed framework for constraining hadronic physics by QCD 
symmetries is chiral perturbation theory (ChPT) |jl|], which provides a systematic expansion 
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in energy for low-energy scattering processes. The degrees of freedom are the Goldstone 
bosons (pions, etc.) and, when appropriate, nucleons. This approach builds in constraints 
due to chiral symmetry without any additional constraints on the dynamics or ad hoc model 
assumptions; physics beyond chiral symmetry is incorporated through constants in the low- 
energy lagrangian, which are usually determined from experiment. Because additional con- 
stants are needed at each stage in the energy expansion, ChPT is predictive only at suffi- 
ciently low energies, where the number of parameters introduced does not overwhelm the 
data to be described. 

The prospects for extending ChPT in a useful way to calculations at finite density are 
unclear at present. On the other hand, the general framework of ChPT has validated the 
principle of resonance dominance of low-energy QCD. In particular, the E"^ coupling con- 
stants in ChPT in the meson sector are well reproduced from a meson resonance lagrangian 
applied at tree level, with the vector mesons playing the leading role 0. Meson dominance 
is also the key principle underlying phenomenological models of nuclei with hadronic degrees 
of freedom, which we consider here. But while the correspondence in the vector channels 
is relatively straightforward because of well-defined resonances, the dynamics in the scalar 
channel is more difficult to identify and to model. 

Within meson-exchange phenomenology, the mid-range attraction between nucleons is 
generally believed to be a dynamical consequence of the strong interactions between two 
pions exchanged with scalar, isoscalar quantum numbers 0. No nearby underlying reso- 
nance at the relevant mass (~ 500 MeV) is evident or, in principle, needed. (Note that in 
ChPT investigations, the scalar resonance is identified with mesons around 1 GeV.) Never- 
theless, this physics is efficiently, conveniently, and adequately represented at the one-meson 
exchange level by the exchange of a light scalar degree of freedom [Q. This light scalar is 
also an essential element of phenomenologically successful mean-field models of nuclei 

These mean-field models are significantly constrained by the bulk properties of finite 
nuclei The question then arises: How should QCD symmetry constraints be manifested 

in these models? There is a long history of attempts to generalize the linear sigma model 
to build models with chiral symmetry; it is almost irresistible to identify the scalar meson 
mediating the mid-range nucleon-nucleon (NN) attraction with the chiral partner of the 
pion. More recently, interest in models realizing the broken scale invariance of QCD has been 
revived. Scale invariance is particularly compelling to consider because of its connection to 
the scalar channel. The breaking of scale invariance by the trace anomaly implies relations 
involving zero-momentum Green's functions of the scalar trace of the energy-momentum 
tensor [see Eqs. (0)-(|^)]; these are called low-energy theorems 0]. If these relations are 



assumed to be saturated by scalar particles at tree level, significant constraints arise on the 
associated scalar potentials (in the chiral limit). We will exploit such constraints in this 
paper. 

In Ref. 0, a broad class of models that attempt to unite successful mean-field phe- 
nomenology with chiral symmetry and the broken scale invariance of QCD were studied. 
Generalizations of the conventional linear sigma model that feature a "Mexican hat" po- 
tential were found to fail generically, even with modifications inspired by the realization 
of broken scale invariance. A significant improvement was found by the Minnesota group 
||TT| when the "Mexican hat" potential is abandoned, and a reasonable description of the 
properties of closed-shell nuclei was obtained. 
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In this paper, we build a different effective model of nuclei by implementing a nonlinear 
realization of cliiral symmetry together with the low-energy theorems of broken scale invari- 



ance. The detailed construction of the full model will be reported elsewhere []I2[ . Our focus 



here is primarily on how vacuum dynamics might be treated in an effective field theory of 
nuclei. 

The role and manifestation of vacuum dynamics is an important issue in any field- 
theoretic description of nuclear matter and finite nuclei. Valence nucleons in the Fermi sea 
interact with each other and also with the QCD vacuum. In turn, the vacuum is modi- 
fied by interactions with valence nucleons. In nonrelativistic models, such effects are never 
dealt with explicitly, but are absorbed implicitly into phenomenological effective interac- 
tions involving only valence nucleons. As a result, the interactions may acquire additional 
density dependence and nonlocalities. In previous relativistic models of nuclear matter in- 
volving a scalar field coupled to the nucleon, vacuum modifications were incorporated in the 
renormalized scalar effective potential [^Q. This in turn affects the density dependence. 

In principle, the one-baryon-loop effective potential contains an infinite number of unde- 
termined coupling constants, which are the coefficients in a polynomial of infinite order in 
the scalar field. In conventional renormalizable models, the nucleon vacuum one-loop correc- 
tion is well defined |]T^ and determines these coefficients, except for the terms of degree four 
and less, which are fixed by a renormalization prescription. However, renormalizable models 
with one-loop corrections do not achieve the phenomenological success of models without 
vacuum terms for the bulk properties of finite nuclei |]14|,|15|. We interpret this failure as a 
phenomenological indication that the vacuum is not treated adequately. 

In previous studies involving nonrenormalizable models, the effective potential is sim- 
ply truncated, usually at degree four, and mean-field theory is applied without considering 
vacuum effects. In this paper, we begin to address the problem of constructing consis- 
tent calculations in effective field theories of nuclear matter and finite nuclei that explicitly 
address the role of the vacuum dynamics. In particular, we show how vacuum loop con- 
tributions are absorbed in the renormalization of coupling constants in the lagrangian in a 
model constrained to satisfy the low-energy theorems of QCD. 

In contrast to the situation in ChPT, we cannot expand in powers of the energy, since 
we are not limited to derivative couplings and light meson masses. We observe, however, 
that the meson fields develop nonzero expectation values (mean fields) at finite density, 
and to begin, we assume that these mean fields dominate the contributions to the energy. 
Successful mean-field phenomenology shows that for densities not much higher than nuclear 
matter equilibrium density, the corresponding mean fields (or nucleon self-energies) are small 
compared to the free nucleon mass (roughly | to | the size); these ratios are therefore useful 
expansion parameters. Moreover, since the derivatives of the mean fields are small for normal 
nuclei, a truncation of the lagrangian at some low order of derivatives is also appropriate. 
(We verify this assertion explicitly later.) The end result is an energy functional for nuclear 
matter and nuclei that contains small powers of the mean fields and their derivatives; nuclear 
phenomenology implies that these fields are an efficient way to incorporate the density 
dependence of nuclear observables. Our objective here is to see how the low-energy behavior 
of QCD constrains the coefficients in this energy functional, particularly with regard to 
contributions from the quantum vacuum. 

The assumption of mean-field dominance also has phenomenological support from Dirac- 
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Brueckner-Hartree-Fock (DBHF) calculations, which indicate that exchange terms and 
short-range correlations do not significantly change the size of the nucleon self-energies 
nor introduce strong momentum dependence (at least for occupied states) Thus 
we have the favorable situation that the mean fields are large enough (compared to nuclear 
energy scales) to dominate the bulk dynamics but small enough (compared to the nucleon 
mass) to provide useful expansion parameters. Going beyond one-loop order systematically 
is an essential issue, but we will leave this as a topic for future study. 

A nonlinear realization of chiral symmetry will be adopted, in which the Goldstone 
bosons (pions) are derivatively coupled to the nucleons. Historically, a linear representation 
(as in the usual linear sigma model) has been favored by model builders, in part because the 
sigma model is renormalizable. In this work, we wish to introduce a light scalar degree of 
freedom, but we do not want to make the restrictive dynamical assumption that this scalar 
is the chiral partner of the pions in a linear representation. By realizing chiral symmetry 
nonlinearly, we are not committed to such assumptions about the scalar degree of freedom. 
In addition, it will be easier to introduce vector mesons in a chirally invariant way that 
manifests the vector-meson dominance of Sakurai [jT8[. Finally, the nonlinear representation 



is more efficient for preserving the consequences of chiral symmetry at finite density when 
making approximations involving pions, because sensitive cancellations are not needed 

As suggested above, the nonderivative terms of the light scalar effective potential can be 
constrained by the low-energy theorems of QCD, so that vacuum effects are "built in." This, 
together with the truncation of our expansion in derivatives and powers of fields, leaves us 
with relatively few parameters, which can be determined by fitting to the properties of finite 
nuclei. 

The paper is organized as follows: In Section II, broken scale invariance is discussed and 
the model is introduced and its renormalization is considered. An approximation scheme 
for nuclear matter and finite nuclei is proposed in Section III and the energy functional is 
derived. Results are given in Section IV. Section V contains some discussion of the results 
and Section VI is a summary. 



II. THE MODEL 

In meson-exchange phenomenology there is a light scalar degree of freedom that simulates 



two-pion-exchange physics in the scalar channel yJSQ] . Here we seek to describe this physics 
by introducing a light scalar field S{x). We do not associate the scalar with a bound state 
or resonance, so we allow S{x) to have anomalous behavior under a scale transformation in 
the effective theory. In particular, when x — >• \~^x, S{x) — >• X'^S{Xx), where d can differ 
from unity and is to be determined phenomenologically. A QCD-inspired scenario that leads 
to such a scalar was proposed by Miransky and collaborators They introduced a light 
scalar generated by dynamical chiral symmetry breaking in QCD, which was consequently 



We are unaware, however, of any proof of the independence of finite-density observables with 
respect to nonlinear field transformations, analogous to the theorem that applies to 5-matrix 
elements [p^ . 
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associated with the quark condensate {qq) and referred to as quarkonium. We will take all 
other fields to have canonical scale dimension. 

While massless QCD is scale invariant at the classical level, this symmetry is broken at 
the quantum level. This breaking is manifested in a nonzero trace of the energy-momentum 
tensor of QCD, which is referred to as the trace anomaly. The QCD trace anomaly in the 
chiral limit is given by 

^/^(a;) ^ -H{x) = {P{g)/2g)G%G'^^'' , a = 1, 2, ■ ■ ■ , iV^ - 1 , (1) 

where G^j, is the gluon field tensor and l3{g) = —{g^/4:8n'^){llNc — 2N{) is the one-loop 
beta function with A^c colors and A^f fiavors. There are remnants of scale invariance, which 
imply low-energy theorems that relate Green's functions involving the trace of the energy- 
momentum tensor H[x) |1T0| : 

i J d^x {0\T[H{x)H{0)]\0) = AH,, , (2) 
z'J d'xd'y {0\T[H{x)H{y)Hm\0) = A^H, , (3) 



where = {0\H\0). 

Effective lagrangians for pure-glue QCD (no quarks) featuring a scalar glueball field x(^) 
("gluonium") that saturates these low-energy theorems at tree level have been considered 
many times 0] . Lattice QCD calculations indicate that the scalar glueball is quite heavy on 
hadronic mass scales, with a mass of roughly 1.6-1.8 GeV [^. Its fate in the real world with 
light quarks is not entirely clear. Here we will generalize the effective gluonium model to 
include the light scalar discussed earlier; this extension was proposed in a different context 
in Ref. pl|. We take the trace anomaly to consist of two contributions, corresponding to 
a vacuum expectation value = H„ + H„. Here H„ is identified with the heavy glueball 



contribution, while i/q is nonzero only when chiral symmetry is dynamically broken in the 
presence of light quarks. One can argue that dominates (which is equal to the gluon 
condensate up to a factor) so that ifq <^ |21[]. How the QCD trace anomaly actually 
separates into the two parts is not explored here, since we will determine ifq by fitting to 
the properties of finite nuclei. Nevertheless, we find that the value of ifq determined in our 
fits satisfies H^^ <^ Hq (see Table |). 

The low-energy theorems involving the trace 0'^{x) of the energy- momentum tensor are 
assumed to be saturated by the scalar gluonium x(a;) and the light scalar S{x). For 
simplicity we adopt a model with no mixing between the scalars. A candidate effective 
lagrangian of the scalars that satisfies the low-energy theorems at tree level in the chiral 
limit is 



23,21 



1 1 



X 



2x 1-d 







d^Sd'^S-Vix,S) 



(4) 



where ai is a real constant, d is the scale dimension of the S{x) field, and the scale-breaking 
potential V is 
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Here xo and 5*0 are the vacuum expectation values of x and S respectively. Notice that ai 
has been introduced so that after expanding the terms in square brackets in Eq. (^), the 
kinetic term for 5* is canonical. The mass of the light scalar 5* is given by ml = 4Hq^/ {d'^ Sq) . 
The scale dimension of the x fi^ld is assumed to be unity. 

One can define the energy-momentum tensor so that the Noether current for scale trans- 
formations is XyO^'^ . The trace of this "improved" energy- momentum tensor correspond- 
ing to the lagrangian in Eqs. (^) and @ is 

dV dV 

o',i^) = sd^+Xg^-mx,s) 

With the dynamics of the scalar field fluctuations governed by the lagrangian in Eqs. @ 
and the preceding trace satisfies the low-energy theorems at the tree level. The 
usual direct demonstration in which the gluonium alone is assumed to saturate the 
Green's functions, involves parametrizing the fluctuations x{^) the exponential form 
X = Xoexp[x(a;)/xo] and substituting this into the low-energy theorems and into the poten- 
tial, V{x, 0) of Eq. (^, to determine vertices. Keeping only tree level diagrams (no loops), 
the theorems then follow. 

To extend the demonstration to the present case, one first notes that the low-energy 
theorems should not depend on how the gluonium field fluctuation x is parametrized. One 
then observes that if x is defined through x = Xo(l — x/XoY, the resulting form for the 
gluonium parts of the trace and the potential in Eqs. (|^) and (|^) become the same as for 
the light scalar, when the fluctuation of the latter is parametrized simply as S = Sq — (p. 
Since there are no couplings between the x and the S fields, the low-energy theorems follow 
directly. 

We can now add to Cg a scale-invariant lagrangian with these scalars coupled to pion, 
nucleon, and vector degrees of freedom. (We neglect pion mass terms at this point.) The 
resulting model would be a candidate model for nuclei that satisfies the low-energy theorems. 
On the other hand, there are many other possible terms allowed, and even Eq. (^ is not 
the most general form involving two scalars. We choose to take advantage of the heaviness 
of the gluonium and the usefulness of an expansion in powers and derivatives of the other 
fields to both simplify and generalize the effective lagrangian. We expect the expansion to 
be valid and useful when applied near normal nuclear matter densities. 

Since the mass scale of the heavy gluonium field (roughly 1.6-1.8 GeV) is significantly 
higher than the scales involved in the nuclear matter problem, the heavy gluonium field 
fluctuations x can be integrated out as in Ref. 0. In particular, we can eliminate x by 
iteratively solving its equation of motion, exploiting the dominance of the mass term over 
powers and derivatives of x- This results in complicated terms involving powers and deriva- 
tives of the other fields, but we can expand these terms. For example, the second term in 
Eq. (H) would become 

^^l + (3,^ + l3,^ + ..] d,(t>d^(t> + Hd,(l^d^4>? + ■ ■ ■ , (7) 

Do Oq J 
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where the /3i are functions of the constants xo, cti, and d. The lagrangian will be much 
simpler if we can truncate this expansion at leading order in derivatives and neglect high 
powers of the meson fields. 

We can follow this prescription to write a general chirally invariant effective lagrangian 
for nuclear matter and nuclei. The ground states of even-even nuclei and nuclear matter will 
be assumed to have good parity, so there is no pion mean field. Thus the pion will not play 
an explicit role in the present discussion of uniform nuclear matter and closed-shell nuclei in 
the Hartree approximation. Nevertheless, we wish to stress the connection to pion physics 
and the underlying constraints of chiral symmetry. Thus, we give an overview of the full 
model in order to motivate the form of the lagrangian and to set the stage for future work. 

We restrict consideration to a low-energy representation of massless, two-flavor QCD. 
The Goldstone pion fields are represented by a chiral phase angle that corresponds to a pure 
chiral rotation of the identity matrix: 

^U = U{x)=exp{i7T{x)-T/U) , (8) 

where C,{x) = exp(?7r(x)-T/2/^), ( a = 1,2, and 3 ) are the Pauli matrices, 7r"(x) are the 
Goldstone pion fields, and /^r = 93 MeV is the pion-decay constant. This parametrization 
and the nucleon representation that follows is conventional; see, for example, Ref. The 
nucleon field is written as 




with p{x) and n{x) being the proton and neutron fields. 

Under chiral transformations of SU{2)i^ (g) 5'f/(2)R, U{x) transforms globally, U{x) — >• 
LU{x)B), where L = exp^iO^ • t) and R = exp{iOFi • t) are x-independent elements of 
SU (2)l and SU (2)r, respectively. In general, the transformation of ^ is local since it depends 
on the pion field p8| , p6| ,^ 



^(x) ^ = mx)h\x) = h{x)^{x)R^ , (10) 

where the second equation defines the SU (2)-valued function h{x) as a nonlinear function 
of L, R, and U{x). Note that h = L when L = R, i.e., in the case of a pure isospin rotation. 
The nucleon field also transforms locally: N{x) h{x)N{x), which implies that nucleons 
mix with pions under chiral transformations. (See Ref. ||2^ for alternative representations 
of the nucleon field.) 

We will incorporate the physics of vector dominance in our lagrangian by introducing 
vector mesons as gauge bosons |jT9|. For simplicity, since we concentrate on the properties 



of nearly symmetric (A^ ~ Z) nuclear matter in this paper, we will not explicitly write down 
the rho and the electromagnetic fields. Thus only the u meson field appears explicitly 



here. We will present a full discussion of the lagrangian elsewhere [|T2[, including how the 
vector mesons are gauged and how vector dominance results. To build chirally invariant 
terms, it is useful to define the vector and axial vector fields 

v,i^) = ~{^^d,^ + ^d,e) , (11) 
a,{x) = -'-{ed,^-^d,e) , (12) 
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which transform as — hv^^h^ — ihd^h^ and ha^j^h^ . The couphng of the pion to the 

nucleon is reahzed through and the covariant derivative 

1^f, = d^ + iv^ + ig^V^ . (13) 

Now we can write the complete chirally invariant lagrangian; all terms not contained in 
Cg are scale invariant. After integrating out x and expanding about Sq, the lagrangian takes 
the form 



C(x) = NiirfV^ - ig^Yli.^ -M + g,4,+ ---)N- jFj^f" 



+ 2 



4 

1 



1 1 / / 1 1 ~ 



where = 1-23 is the axial coupling constant, gs {g^) is the light scalar (vector u) coupling 
to the nucleon, the u field strength tensor is F^^ = d^V^ — d^V^, and rj and C are real 
constants. 

Several features of this lagrangian are of interest: 

• We have combined terms after expanding and have rewritten the coefficients, where 
appropriate, in terms of physical masses. Note that the nucleon mass M has contri- 
butions from the vacuum expectation values of both scalars; we do not assume that it 
comes entirely from the light scalar (although this possibility is not excluded). 



The combination of the uj mass term and the pion kinetic term in Eq. ([Ij) appears 



naturally, if we assume the vector mesons to be gauge bosons [19,12 



• The original separation of the lagrangian into a scale-invariant piece and a scale- 
breaking piece, in which the latter involved only the scalar fields, is now largely hidden 
because the x dependence is not exphcit and we have expanded about 5*0. Nevertheless, 
there is a remnant for our purposes here: the scale-breaking potential of the light scalar 
[the last term in Eq. (0)], which is not changed by the elimination of x- (Recall that 
the X and S do not mix.) Thus the low-energy theorems still protect the form of 
this potential, which places constraints on vacuum loop renormalizations, as discussed 
below. 

• We have omitted many higher-order terms, as indicated by the ellipses, which represent 
higher powers of fields and their derivatives. Only Yukawa couplings to the nucleon 
fields are kept, based on the phenomeno logical dominance of one- meson exchange and 
the implicit elimination of heavier fields. (So NNcfP' terms, etc., are omitted.) Higher- 
order terms with meson fields should give numerically small contributions (in nuclei) 
or can be absorbed into slight adjustments of the other parameters. Some explicit 
justification for these claims is given in the results below. 

The lagrangian in Eq. (|1^ is written with renormalized coefficients. Counterterms are 
not written explicitly, but are implied. In particular, these counterterms include all powers 
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of the scalar field, not just terms up to 0(0^), as in a renormalizable model. To understand 
how these counterterms are fixed, we start by integrating out the baryon fields at zero density 
and temperature. The result is a fermion determinant that contributes to the meson action 
as an additive term given by 

Sm[<P,V^] = J d^xAd = -tTT\nK{0) , (15) 

where "Tr" indicates a trace over spacetime, spin, and isospin, and the kernel K{fi) is defined 
in coordinate space by 

{x\K{fi)\y) = [z^d, - g^rV^ix) + fijo - M + g,(P{x)]S\^ - v) ■ (16) 

The introduction of the chemical potential /i is for later convenience, and baryon countert- 
erms, which are needed beyond one-loop, are suppressed. Note that no approximation has 
been made at this point; 5'fd is a functional of the dynamical fields and that still must 
be integrated over in a path integral, for example. The techniques for expanding a determi- 
nant in powers of derivatives can be found in Ref. |2^; see also the heat-kernel method in 
Ref. p6|. The expansion of Eq. (|15D in a renormalizable model has been discussed in Ref. 



We first focus on the nonderivative terms, which can be obtained from Eq. (|T5D by 
treating the fields as constants and by expanding the logarithm in a power series in the 
fields. Baryon number conservation implies that for the vector field, only its derivatives can 
appear in the expansion. Thus the nonderivative part of £fd is an infinite polynomial in 0; 
for example, at the one-loop level, 

AdM = z tr In G°(A:) + . E ^[^7s0(x)]" tr [G\k)T . (17) 

Here we have regularized dimensionally to maintain Lorentz covariance and baryon number 
conservation, "tr" denotes a trace over spin and isospin only, and 

is the free baryon propagator. Beyond one loop there are additional terms in the coefficients, 
including baryon counterterm contributions. 

The polynomial in of Eq. (^) must be combined with the corresponding counterterms; 
in this way the vacuum contributions are absorbed into the renormalization of the scalar 
polynomial. If one insists that the low-energy theorems be satisfied at tree level in the meson 



fields, the end result for the scalar potential should be of the form in Eq. (T|), where the 
couplings are renormalized. (Note that this potential can be expanded as a polynomial in 0, 
with all coefficients determined by ifg, 5*0, and d.) One never has to explicitly calculate any 
counterterms or evaluate Eq. (0); when we write down the scalar potential, the nucleon-loop 
effects have already been taken into account. Furthermore, although we have illustrated the 
renormalization by evaluating nucleon loops only, any additional baryonic degrees of freedom 
in the lagrangian would be treated analogously and the final result will be the same. Thus 
the phenomenological fitting of parameters accommodates a general characterization of the 
vacuum response. 
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The renormalization of the derivative terms is analogous except that we do not have 
low-energy theorems to reduce the number of renormalized coupling constants. We note, 
however, that each additional derivative is accompanied by an inverse power of a typical 
scale in the problem, which is the nucleon mass here. Experience with mean-field mod- 
els of nuclei also suggest that the derivatives of the mean fields are small (for example, 
I V0/0I ^ 100 MeV). Thus if we assume mean-field dominance, such that fluctuations around 
the mean fields are small, and the naturalness of the coefficients in the derivative expansion 
(see the discussion in Section V), we can truncate the derivative terms at some tractable 
order. In this work, we will stop at the lowest order for the derivatives. Thus we have 
only a few unknown renormalized constants (parameters), which are determined by fitting 
to experiment; in our case, we will use finite-density observables. 

At finite density, we work in the grand canonical ensemble through the introduction of 



a chemical potential fx I^Tf . We consider only zero temperature in this work, which allows a 



simplified discussion. The relevant lagrangian density is now 

C'{x, /i) = C{x) + fxNjoN . (19) 

Here the effective action of C is associated with the thermodynamic potential fl of the 
system, instead of the energy. The energy follows from 

E = n + fxB , (20) 

where 

—I (-) 

is the baryon number of the system. 

Now we integrate out the baryon field as at zero density. The result is the fermion 
determinant at finite density (or chemical potential), —i Tr lni^(/i), to which we can add and 
subtract the fermion determinant at /i = 0, — zTrlnK(O).0 The added term — iTrlnii'(O) 
combines with the counterterms exactly as described above so that the renormalization goes 
through as before. Note that it contains the same dynamical scalar and vector fields as the 
fermion determinant at /i. The remaining combination 

-iTi\nK{fi)+iTi\nK{0) (22) 

is an explicitly density-dependent piece (it vanishes for fi = 0), which is finite if baryon 
counterterms are included in K{fi). (This combination is evaluated in the Hartree approxi- 
mation in the next section, for which the baryon counterterms are not needed.) Once again 
the scalar potential in the form shown in Eq. (0) is left intact; the only difference is that 
the scalar field now acquires a different expectation value due to the presence of valence 
nucleons at finite density. 



^We assume that /i = still separates the positive-energy levels from the Dirac sea. This will be 
the case if the density is not too high. 
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III. FINITE NUCLEI AND NUCLEAR MATTER 



To perform a realistic calculation, we need a good starting approximation. Since our 
focus here is on bulk nuclear properties and on single-particle spectra, we assume that the 
mean meson fields dominate the dynamics, and we expand the finite-density thermodynamic 
potential around the mean fields. The lowest-order result (Hartree approximation) is ob- 
tained by replacing all the meson fields by their mean values, and this will be the starting 
point of any systematic approximation for treating the fluctuations. 

The thermodynamic potential for nuclei in the Hartree approximation is given by 

y dxo n = 2 Tr In Z(/i) - iTT\nK{0) ~ J d'^x (x) , (23) 

where the baryon kernel in coordinate space is now 

{x\K{fi)\y) = Jo[^^o + - h{^)]6^^\x - y) . (24) 

The single particle hamiltonian h is 

/i(x) = -2a • V + g^Voi^) + [3{M - gM^)) , (25) 

with /3 = 7o and ct = 707, and the static scalar and vector mean fields are denoted by </'o(x) 
and Vo(x). The contribution from the meson fields is 



+ 2( VVy^ + - I 1 + m%' + jiCig-VoY . (26) 

Note that K{fi) is diagonal in the single-particle basis ijja{'x)e^'^^° , where V'a(x) are the 
normalized eigenfunctions of the Dirac equation with eigenvalues P^,|5[|: 

hM^) = ^a^a(x) , J d^xV-t (x)^,(x) = 1 . (27) 

From a path integral formulation, one can see that the appropriate boundary condition or 
ie prescription for evaluating the baryon kernel is — (1 + ie)uj. 
From Eq. (|23|) one can now obtain, after a Wick rotation, 

= - V / — [ln(-iw + fx- E^) - ln{-iuj - E^)] - [ d^x 

= -Y,{fi-E^)[e{fi-E^)-ei-E^)]- f d^xu^ 

a •' 
occ „ 

= -Y.{li-E^)- j d^'xU^ . (28) 



Here we have used 



E^(-^") = E^(^a) = E^ (29) 

r\i ^ 
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which is vahd when = separates the nucleon levels from the antinucleon levels. The 
summation superscript "occ" means that the sum runs only over occupied states in the 
Fermi sea. Moreover, using Eqs. ( pO]) and (0), we find 

occ 

b=y: [0{f^ - - ei-E^)] = 5: 1 , (30) 

a a 
occ „ 

E = J2E^- d'xU^. (31) 

a ■' 

We emphasize that the final sum over only occupied (valence) states is not the result of a 
vacuum subtraction, as the term with n = still contains the background fields, which must 
be determined self-consistently. The true vacuum subtraction was performed earlier when 
we derived the renormalized f/m- 

The equations for the mean fields are obtained from extremizing the energy functional 
with respect to 0o(x) and Vo(x). From Eqs. ( p5[ ) and (|27| ) one finds 

(5K S 



6o(x) 500 (x) 



Oh 6 f 

= (x)— 7/;„(x) + E^-—^ / d'y£iy)My) 
O(po c'0o(x) J 

dh 

= V^i(x)^V'«(x), (32) 
and a similar expression for the variation with respect to Vq; evaluating the derivatives yields 

^ occ occ 

J2Ea = -^s$^^a(x)?/'a(x) , (33) 



^OlXy „ a 

^ occ occ 



J2E^=g^Y.^i{^)M^) ■ (34) 



5^o(x) ^ 

Upon applying these results to Eq. (|3lD , one obtains the mean-field equations: 

occ 

- V^0o + ml(j)o = gsJ2 ^a(x)^a(x) 

a 

+ m^(l)o + m^Soil - — ] In 1 - — + — m^\/o , (35) 



5*0/ V 5*0/ 25*1 



-VVo + mlVo = ^.^V^l(x)V^„(x) 



- v'^rnlV, - hgX' ■ (36) 

Note that we have added an explicit mass term to each side of the scalar field equation to 
put it in a form that can be solved with conventional numerical techniques [^|. The vector 
mean-field equation is actually a constraint since the time component of the vector field is 
not a dynamical degree of freedom. (See below for further comments in the case of nuclear 
matter.) 
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This lowest-order result (Hartree approximation) is similar to that obtained from conven- 
tional derivations of relativistic mean-field models in which one-loop vacuum corrections are 
simply neglected. We emphasize, however, that we are not merely presenting another mean- 
field model; the vacuum effects are incorporated and systematic improvement is possible 
(in principle). Rather, we consider our procedure a justification for the phenomenologically 
successful mean-field approach. 

The energy density for uniform nuclear matter in the Hartree approximation can be 
obtained from the preceding results by observing that the single-particle energy eigenvalue 
becomes 



E{k) = g^Vo + + M*' , (37) 

where M* = M — gs4>o, and 0o and Vq are now constant mean fields. The energy density S 
becomes 



1, 1 



1/. , 0O\ 2T.2 , 7 f''^ 



It 1 



f [M-, p.] = -miSld'-{{l - f 1„(1 - - i] + i} + g^.pM - ^((SM) 



i(l + r,|) ,nX' + p2^/ 'Sk Vk' + M'^ , (38) 

where kp is the Fermi momentum defined by = gvVo + \/kp + M*^, and pa = 7A;p/(67r^) 
is the baryon density. The spin-isospin degeneracy 7 = 4 for nuclear matter and 7 = 2 for 
neutron matter. 

The equation that determines Vq can be obtained either from the Euler-Lagrange equa- 



tions or by using Dirac's procedure |23[, with the result 



= (1 + V^) mlVo + l(gX ■ (39) 

This equation can also be obtained from Eq. (|38|) by setting {dS / dVo) p^^^M* = . Note, 
however, that this is not a minimization condition for S. In fact, the Vq obtained from 
Eq. ( |39D corresponds to a local maximum of the energy density. Equation (pOf), like Eq. (p6[), 
is a constraint equation for Vq, which is not a dynamical variable. 

The energy density at a given baryon density is found by using Eq. (^) to eliminate Vq 
from S in Eq. (|38D and then by minimizing the resulting finite-density effective potential 
with respect to M*. The effective potential at fixed baryon density is shown in Fig. 
Notice that in contrast to the conventional one-loop approximation (relativistic Hartree 
approximation or RHA [^,^) in renormalizable models, the finite-density effective potential 
of our truncated model is meaningful only when \gs4>o\ is sufficiently small that higher-order 
terms can be neglected. Similar considerations apply to the solutions of Eq. (|39D . (See 
Section V for further discussion.) Parameters can be chosen so that nuclear matter exhibits 
saturation at the empirical point; one approach to determining the parameters is discussed 
in the next section. 



IV. RESULTS 

To test the utility of the model, we must see if it can successfully describe finite nuclei 
0. The basic features we seek to reproduce are the nuclear charge densities (including 
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FIG. 1. Finite-density effective potential 8 from Eq. (p8|), plotted as a function of M* (solid 
line). Vq is eliminated for each M* using Eq. (^). Parameter set Tl is used and = 1.30 fm~^. 
Results for other parameter sets and other densities are qualitatively similar. Also shown are the 
analogous potentials for the Walecka model RHA (dotted line) and the nonlinear parameter set 
B from Ref. 0] (dashed line). 



the observed flatness in heavy nuclei), the characteristics of the single-particle spectrum, 
and the bulk binding-energy systematics. Relativistic mean-field models unconstrained by 
QCD symmetries have been successful in reproducing these properties for nuclei across the 
periodic table. 

The Hartree equations for finite nuclei in our model were given in Section III, but only 
isoscalar mesons were discussed. To make realistic comparisons to experiment, we must 
include the p and the Coulomb interactions. Here we simply introduce the p and the photon 
as in Ref. , except that we also include a coupling between the p and the scalar 0, exactly 
as for the u [see Eq. (^81)]. A more complete treatment of the isovector mesons will be 
presented elsewhere |T^ . 

We take the nucleon, cj, and p masses as given by their experimental values: M = 
939 MeV, rriy = 783 MeV, and rrip = 770 MeV. We then fit the rest of the parameters {qs, 
(7v, dp, Ci So, and d) to the binding energies, the charge radii, and the spin-orbit 
splittings of the least-bound proton and neutron in ^^O, '^''Ca, and ^'^^Pb, as well as to the 
charge density of ^^O at r ^ Ifm. An optimization process similar to that of Ref. p5| is 



used. Here we are principally interested in showing that a good fit to properties of finite 
nuclei can be achieved; Table | lists three such parameter sets (Tl, T2, and T3). In set Tl, 
d is an optimization parameter, while it is fixed (arbitrarily) in sets T2 and T3 to illustrate 
the range of possible d. In a future paper, we will study in more detail the regions of the 
parameter space that produce a reasonable fit and examine which conditions are important 
in determining individual parameters. 

We have calculated ^^O, ^°Ca, and ^ospb for these parameter sets and for a representative 
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TABLE I. Parameter sets from fits to finite nuclei. The vector masses are mv = 783 MeV and 
nip = 770 MeV; the nucleon mass is M = 939 MeV. Values for Sq, the scalar mass mg, and hI/'^ 
are in MeV. Note that ml = AHj{(fSl). 



Set 


9l 


TTls 


9l 




5*0 


c 


V 


d 


0-1/4 


Tl 


99.3 


509. 


154.5 


70.2 


90.6 


0.0402 


-0.496 


2.70 


250. 


T2 


96.3 


529. 


138.0 


69.6 


95.6 


0.0342 


-0.701 


2.20 


236. 


T3 


109.5 


508. 


178.6 


67.2 


89.8 


0.0346 


-0.160 


3.50 


283. 



TABLE IL Binding-energy systematics for the model proposed here (sets Tl, T2, and T3), for 
model B from Ref. Q, and for the point-coupling (PC) model of Ref. ||3^. Binding energies per 
nucleon are given in MeV. 



Model iOQ^ 208p^3 

^fl 7^99 8^61 7.91 

T2 7.94 8.55 7.89 

T3 7.95 8.53 7.91 

B 7.82 8.35 7.62 

PC 7.97 8.58 7.87 

exp't 7.98 8.55 7.87 



mean-field model (set B from Ref. 0). Bulk binding-energy systematics are summarized in 
Table ^ and rms charge radii are summarized in Table |ITT|. For comparison, we also include 



results from the point-coupling model of Ref. |35|. The binding energies include center-of- 



mass corrections as in Ref. 0. We show charge densities and single-particle levels for ^osp]-, 
in Figs. and ^, and charge densities for ^^O and ^°Ca in Figs. ^ and ^. The charge densities 



are determined from point-proton densities following the conventional procedure [32|, which 



folds them with a phenomenological proton form factor. Form factors generated within the 
model itself, originating from vector dominance physics, will be considered elsewhere. 

The fits to nuclear charge radii, binding energies, and spin-orbit splittings are quite good. 
The only deficiencies in the sets illustrated here are some small deviations from experiment 
in the charge densities. Changes in the optimization procedure can improve the agreement 
of the charge densities at the cost of worsening slightly the agreement with empirical binding 



TABLE III. Rms charge radii (in fm) for the model proposed here (sets Tl, T2, and T3), for 
model B from Ref. Q, and for the point-coupling (PC) model of Ref. 



Model iOQ^ 208p^j 

^Tl 2?73 3^7 5.56 

T2 2.72 3.47 5.56 

T3 2.72 3.48 5.57 

B 2.74 3.48 5.56 

PC 2.73 3.45 5.51 

exp't 2.74 3.47 5.50 



15 
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FIG. 2. Charge density of ^'^^Pb. The sohd hne is taken from experiment [37|. Charge densities 
are shown for a successful mean-field model (model B from Ref. |§) and for the three parameter 
sets from Table I. 



energies. 

A good reproduction of the spin-orbit force in finite nuclei necessarily leads to large scalar 
and vector mean fields in the interiors of the nuclei or in nuclear matter. In particular, as 
discussed many times (recently by Bodmer P]), vector and scalar fields of roughly 250- 
300 MeV are needed to reproduce the observed spin-orbit splittings in the least-bound levels 
(and also the deformations in light, axially symmetric nuclei [0). While these fields are 
large on the scale of the nuclear binding energies, \g^VQ\/M and \gs(t>o\/M and their gradients 
in finite nuclei are relatively small; thus, these remain useful expansion parameters. This 
justifies our truncation of the energy density at small powers of the meson fields. While 
it is possible in principle to add additional monomials in the fields (with undetermined 
parameters), the quality of the present fit makes it unlikely that there is much to be gained 
by this. 

The scale dimension d of the light scalar field was found to be about 2.7 when d was 
included in the optimization. Note that the canonical dimension would have d = 1. Changes 
in the optimization procedure or a relaxation in the goals of the fit allow for a considerable 
range in d (sets T2 and T3 are examples), but it does not seem possible to find a reasonable 
parameter set with d < 2. Thus the introduction of an anomalous dimension for the light 
scalar degree of freedom is an essential feature for the phenomenological success of our model. 

Experience with a broad class of relativistic mean-field models shows that models that 
successfully reproduce bulk properties of finite nuclei share characteristic properties in in- 
finite nuclear matter |p. These properties are the equilibrium binding energy and density, 
the compressibility K., and the value of M* jM at equilibrium. One further condition, that 
the light scalar mass ~ 500 MeV, is needed to ensure reasonably smooth charge densities 
and good surface-energy systematics. If we calculate nuclear matter with the parameter sets 
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EXP'T 



FIG. 3. Predicted proton single-particle spectra for ^'^^Pb using the parameter sets from Table I. 
Only the least-bound major shell is shown. The leftmost values are from experiment, model B is a 
successful mean-field model from Ref. Q, and model PC is the point-coupling model of Ref. |35|. 
Note that the I/19/2 level is an unoccupied state. 
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TABLE IV. Nuclear matter saturation properties for the model proposed here (sets Tl, T2, 
and T3), for model B from Ref. and for the point-coupling (PC) model of Ref. |3^. Values 
are given for the binding energy per nucleon (in MeV), the Fermi momentum kp (in fm^^), the 
compressibility K (in MeV), the bulk symmetry energy coefficient 04 (in MeV), M* /M, and gyVo 
(in MeV) at equilibrium. 



Model 


E/B - M 




K 


04 


M*/M 


gvVo 


Tl 


16.2 


1.30 


194. 


39. 


0.60 


302. 


T2 


16.3 


1.29 


240. 


40. 


0.61 


298. 


T3 


16.1 


1.29 


244. 


34. 


0.61 


297. 


B 


15.8 


1.30 


220. 


35. 


0.63 


277. 


PC 


16.1 


1.30 


264. 




0.58 


322. 



in Table |l[ we find good agreement with values found in investigations with unconstrained 
mean-field models (see Table [TVI ). In particular, the saturation density corresponds to a 
Fermi momentum of about 1.3 fm~^, and the binding energy per nucleon at saturation is 
about 16 MeV. The compressibility is less well determined (190-250 MeV). The nucleon ef- 
fective mass M*/M ^ 0.60 and the scalar mass mg is just over 500 MeV. We emphasize that 
these values are obtained after a fit to finite nuclei only. 



V. DISCUSSION 

We can relate the phenomenological success of the model proposed here to the char- 
acteristics of successful relativistic mean-field models of finite nuclei. A key feature is the 
logarithmic potential for the scalar field, which allows for relatively weak nonlinearities and 
the dominance of the cubic and quartic scalar terms, with the values of the scaling dimen- 
sion d used here. In contrast, chiral models with a Mexican hat potential have large cubic 
and quartic terms, which preclude a good fit to bulk nuclear properties Bodmer p has 
shown that nuclear matter properties that lead to good predictions for finite nuclei can be 
achieved if one adds to (small) cubic and quartic scalar terms a term that is quartic in the 
vector field (here with coupling (). Thus our model has all of the ingredients needed to 
allow a good fit through optimization. In addition, adjustments can be made through the 
scalar- vector coupling 77. 

Note that the scalar-vector coupling and the quartic vector self-coupling can be used to 
define an effective, density-dependent mass of the vector meson at the mean-field level. For 
example, one can use the second derivative of the lagrangian with respect to the vector field. 
For the model parameters in Table |, the two contributions largely cancel, so that the vector 
effective mass m* is essentially independent of density. This is in contrast to the universal 
scaling hypothesis of Brown and Rho [^, which predicts m*/mv = M*/M. 

We have excluded many terms from our model: higher-order polynomials in the vector 
fields and mixed scalar- vector terms, non- Yukawa couplings to the nucleon, derivative terms, 
and so on. In retrospect, were we justified in neglecting them? An analysis of mean-field 
models |§,|12| implies that one can identify dimensionless ratios that can be used to set the 



scale of individual contributions to the energy. For example, one can rewrite the scaled 
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energy density of nuclear matter, £/M'^, in terms of the dimensionless ratios q^Vq/M and 
gs4>Q/M, which then become our finite-density expansion parameters. 

Moreover, an important assumption in applying effective field theories, such as chiral 
perturbation theory, is that the coefficients of terms in the lagrangian are "natural," i.e., of 
order unity, when written in appropriate dimensionless units. This assumption makes the or- 
ganization of terms through a power-counting scheme useful, because one can systematically 
truncate the expansion when working to a desired accuracy. We have proposed an analogous 
concept of naturalness for the finite-density problem, which will justify the neglect of higher 
derivatives and powers of the fields when applying Eq. (|I^) to nuclei. For example, if one 
expresses the nuclear matter energy density in terms of the scaled field variables written 
above, one finds that the ratios 



c 



(40) 



should all be of roughly equal size for our expansion to be "natural." One can verify that 
the values in sets Tl, T2, and T3 satisfy this condition. 

To examine the size of the scalar self-interactions, one expands the logarithmic potential 
in Eq. (p^) with the result 



1 



2\glM^ 



1 



1 



(3d - 8)m2 



?>\Mgl dSo 



1 (llrf2 _48rf + 48)m,2 



+ 



+ 



M 2(25rf3- l40c/2 + 240d- 128; 



5\gl 



(dSo) 



r M2 2(137d^ - 900rf3 + 2040^2 - 1920rf + 640)m2 



(dSo) 



(41) 



where $ = gs4>o/M. The coefficients in square brackets give the combinations that should be 
compared to those in Eq. (^OD , and one can verify that these are also natural for parameter 
sets Tl, T2, and T3. It is interesting that for set Tl, in which d is an optimized parameter, 
the scaled coefficients are extremely small due to nearly complete cancellations among terms 
in the polynomials in d. (Further discussion of these issues is given in Ref. |I^.) 

Further support for the naturalness assumption comes from extending the model to 
include (p'^V^V'^ and {V^V^Y terms and then repeating the optimization. The new fit is 
very close to the fit obtained without these terms. Furthermore, contributions to the energy 
from the new terms are less than 10% of those from the old terms at nuclear matter density, 
and the old coefficients change only slightly in the new fit |jl2[. Thus contributions from the 
higher-order terms can be absorbed into slight adjustments of the coefficients in Eq. (0). 

The astute reader will note that if our naturalness assumption is justified, we could 
construct a variation of our model without the constraints of the low-energy theorems of 
broken scale invariance. Indeed, at nuclear matter density, numerics alone would let us 
truncate the scalar potential, and the same arguments about renormalization apply, so that 
vacuum effects are still built in. This explains the success of previous relativistic mean-field 
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models of nuclear structure and illustrates the power of the assumption of naturalness. Here 
we note that the scalar potential constrained by the low-energy theorems actually provides 
some justification for naturalness, which we simply assume is valid for higher-order and 
derivative terms. Thus one should be cautious in drawing strong conclusions about the role 
of broken scale invariance when applied in models that are restricted to moderate nuclear 
density. 

How widely can our model be applied? A prime motivation for developing relativistic 
models of nuclei and nuclear matter is to extrapolate to extremes of density and temperature 
1^. Such conditions can be reached experimentally in relativistic heavy- ion collisions. One 
hopes that the calibration of such models to observables at ordinary nuclear densities and 
zero temperature, in conjunction with constraints from QCD symmetries, will permit reliable 
extrapolations. 

Unfortunately, our framework of mean-field dominance, naturalness, and the truncation 
at small powers of the fields and their derivatives, which limits the number of parameters at 
ordinary nuclear densities, is bound to break down as the density increases. With increasing 
density, we will find increasing mean fields and expansion parameters that are no longer 
small. Thus we become increasingly less justified in ignoring the effects of higher-order 
terms, and the calibration at nuclear matter density becomes less and less of a constraint. 
The limits of reliable extrapolation are not clear, but one should certainly be cautious in 
applying models like ours much above nuclear matter density. Nevertheless, the utility 
of an accurate relativistic mean-field model for nuclear structure and reactions, which is 
compatible with the low-energy behavior of QCD, should be obvious. 

We close our discussion with some interesting observations. From Table |, one sees that 
5*0, the vacuum expectation value of the light scalar field S, is close to the experimental value 
of (93MeV). Furthermore, the scalar coupling constant gs is close to g-^/gA- If we forget 
for the moment complications from requiring terms to be scale invariant, it is tempting to 
say that the model has a preference for the nucleon mass to be generated entirely from 
the vacuum expectation value of S. That is, if the only nucleon coupling to scalar fields is 
ggNNS, then we recover the empirical nucleon mass and the Goldberger-Treiman relation 
from the fit values of the other parameters. This scenario is also consistent with Miransky's 
model, in which the light scalar (quarkonium) is associated with the quark condensate, 
and with QCD sum rules, which associate the nucleon mass predominantly with the quark 
condensate. It is premature to do more than to point out these results, but the coincidence 
of numbers certainly merits further investigation. 

VI. SUMMARY 

In summary, we have introduced a new model for nuclear matter and finite nuclei that 
realizes QCD symmetries at the hadronic level. In particular, the model incorporates chi- 
ral symmetry, broken scale invariance, and the phenomenology of vector dominance. An 
important feature is the light scalar degree of freedom, which is given an anomalous scale 
dimension. The renormalized scalar potential is constrained by the low-energy theorems of 
broken scale invariance. Vacuum loop effects are absorbed into the renormalized parameters, 
which are determined by fits to hadron masses and finite-density observables. 
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The truncation of the model lagrangian is based on mean-field dominance and the identi- 
fication of expansion parameters that are reasonably small at nuclear matter densities. Due 
to the characteristics of the constrained scalar potential, we adopt a "naturalness" assump- 
tion, which justifies the truncation. The parameters of the truncated model are identified 
by an optimization procedure designed to reproduce bulk properties of finite nuclei. Good 
fits are obtained, which also lead to very reasonable nuclear matter properties. The scale 
dimension of the scalar field comes out greater than two, but is not tightly constrained by 
the fit. 

It is important to emphasize what we have learned about the relationship between effec- 
tive (hadronic) theories of QCD and successful relativistic mean-field phenomenology. The 
vacuum dynamics of QCD is constrained by the trace anomaly and the consequent low- 
energy theorems of QCD. At the level of hadronic fields, this physics manifests itself in the 
scalar-isoscalar sector of the theory. We have proposed that this sector can be divided into a 
low-mass part that is adequately described by a scalar meson with anomalous dimension and 
a high-mass part that is "integrated out," leading to various couplings among the remaining 
fields. We believe this latter characterization is quite general and independent of the details 
of the high-mass part of the scalar sector. Nevertheless, whereas the realization of the Gold- 
stone boson dynamics is well known {i.e., chiral perturbation theory), as is the dynamics of 
the vector sector {i.e., vector- meson dominance), little is known about the precise form and 
magnitudes of the nonlinear couplings originating from the scalar degrees of freedom. We 
find that our primary source of information on this dynamics comes from nuclear structure 
physics, which provides strong constraints on this sector of the theory. 

In subsequent work, we will further explore the parameter space that leads to good fits to 
nuclear properties and identify the observables that constrain individual terms. We will also 
investigate the chiral properties of the model and study the implications of vector dominance 
for nuclear observables. Work to extend the model beyond the one-baryon-loop level in a 
manner consistent with conservation laws and Ward identities is in progress. Finally, we 
will continue the development of the naturalness concept for finite-density systems. 
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